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Sound understanding of the moling efficiency of supercritical hydrogen is crucial to the 
thvd o pa sa t «f high pressure thrust dumber* ter regenerotivety ceded LOX/WI* redid 
engines. This paper examines historical heat transfer correlations for supercritical hydrogen 
wid the effect* ef eecertetetfe* hi hydrogen property date. It h shewn that uncertainty due 
to property data alone can he as high as 10%. Previous heated tube experiments with 
supercritical hydrogen are summarised, and data teem a number ef heated tube 
experiments are analysed to evaluate conditions for which the available correlations are 
vattd. 


Nomenclature 

Cp ■ constant pressure specific best 

D * inner diameter of tube 

h * convection coefficient 

k ■ thermal conductivity 

Nm s Nusselt number 

Pr * Prandtl number 

Pr ~ reduced pressure 

Re ■ Reynolds number 

T « temperature 

Tr * reduced temperature 

V * uncertainty 

V ■ velocity 

x » distance downstream from entrance 

x * any property 

M * dynamic viscosity 

v * kinematic viscosity 

p « density 

a “ standard deviation 

Subscripts 

b ■ property evaluated at bulk temperature reference 

c * critical condition 

eak ■ calculated value 

* experimental value 

/ ■ property evaluated at film temperature reference 

rqt * reference temperature type 

j * condition at Inner-wall surface 

0 4 » property evaluated at 04 temperature 

f 91 property evaluated by integral method 
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I. Introduction 

A nalytical performance prediction* and design margin* for regencrutively corded liquid rocket engines utilize 
correlation* to predict die heat transfer from die combusting gas to the thrust (dumber wall and from the wall 
into toe cooling Quid. These semi^mpirtoat heat transfer correlations are highly dependant on toe fluid transport 
properties and the conditions at which the fluid is flowing. Although a number of coolant-side correlations have 
been proposed for hydrogen acting as toe coolant, toe ability of toe correlations to accurately predict toe heat 
transfer over a wide range of operating conditions ha* typically been limited. 

In addition to toe transport properties, heat transfer is highly dependant on cooling channel geometry Efforts are 
currently being made to produce accurate correction factors that can be applied to existing heat transfer correlations 
to account for toe effects of curvature, roughness, asymmetric heating, entrance effects, and Ugh aspect ratio 
cooling channels for improved thrust chamber design. For these correction factors and new correlations to have any 
merit, it is important to understand toe limitations of these correlations for toe simplest case: a straight, uniformly 
heated, circular tube. This paper examines existing correlations for straight tube geometries using supercritical 
hydrogen as the coolant. The uncertainty on toe heat transfer correlations due to inherent uncertainties in toe 
equation-of-state and transport properties of supercritical hydrogen are evaluated, Previous experimental data ia also 
analysed with existing published correlations over a wide range in toe supercritical region to determine toe regions 
of validity for each correlation, 


H, Kent Transfer Correlations and Studies 


In typical rochet applications, hydrogen enter* toe coolant 
channels as a liquid above toe critical pressure (1.3 MPa) and 
at temperatures as low as 23 K, which is below toe critical 
temperature of 33 K High-pressure thrust chambers using 
hydrogen as a regenerative coolant operate at very high (>10*) 
Reynolds number flow w w ltti^ is in the passages-* As 

the hydrogen continues through toe channels, it is quickly 
heated to above the critical temperature, where it becomes a 
supercritical gas. as shown in Figure 1. The heat transfer 
mechanisms occurring in these regions are supercritical bquid 
and supercritical gas forced convection*. 

Numerous expe rim e n tal studies have been performed for 
hydrogen flowing supereritieatty in uniformly heated circular 
tubes. The flow conditions and geometries for many of these 
studies we summarized in Table 1 From these studies, most of 
which took place in toe 1960s. came a number of correlations 
to describe the heat transfer characteristics of supercritical 
hydrogen in cooling tubes. These correlations are all 
principally based on the conventional Nusselt-type equation of 
Dittus-Boeher, where the Nusselt number is a function of toe 
Reynolds and Prandtl numbers: 1 



Figure 1: Reduced Pressure-Temperature 
Diagram 


Mr » 0.023 Rc‘*Pr M (1) 

In I960, McCarthy and Wolf proposed a modified version of this equation, 4 

Mt « Q-02SRV* P (2) 

where the subscript b denotes that tire properties efevaluated at the bulk temperature of the fluid and T,jT M ia toe 
mner-waU-to-bulk temperature ratio. Hendricks, Summeau and Friedman in 196$ replaced the bulk condition and 
proposed toe correlation, 1 

Nuj • 0.021 Re/'Pr/ 4 (3) 
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where the properties are evaluated at the film t emperatu re. 


h 


Z^ZL 

2 


(4) 


Hess and Kim dan in 1 965, proposed a modified version of the film temperature-type correlation that included 
wall-to-bulk viscosity ratio, 4 


Nu f »0,Q208Re/ , Pr/ J, (l+Q,Q1457v' J /v*) (5) 

Miller, Sender and Ttebes (1965) proposed using a new reference temperature, and developed die correlation, 7 


Nu m * Q,Q2Q8(Re M f (Pr 94 ) M (1 * 0,00983 r, /v * ) 


( 6 ) 


where 


r M -r ( +o,4(r ( -r t ) 


(7) 


This approximately constant reference temperature is equivalent to Deissler’s comprehensive theory which describes 
forced convective heat tranifer for turbulent flow at supercritical gas condition* sufficiently above the critical 
temperature. The use of hulk temperature physical properties with a T t fT k eorreetioo factor, as in the McCarthy & 
Wolf correlation, has also been shown to be equivalent to the Deissler theory, 3 

Taylor, working in 1968*1970 in support of the NHRVA and Fhoebus-2 nuclear rocket engines, added a 
eorreetioo factor to the McCarthy & Wolf correlation to improve the heat fransfer along the length of the tube >l? 

* 0,025R<* (8) 

In 1 973, Sehact and Quentmeyer proposed a new type of reference condition 15 such that, 

Mr, * 0.025 Re^Fr 44 (9) 

This correlation uses an Integral reference type, Property x is evaluated such that, 

x^^VfxtF)^ (,0) 

’• h n 

These correlations have traditionally been used to calculate die heat transfer to supercritical hydrogen cooling 
tubes, and are all based on the Dittus>Boelter4ype equation for Mr , Ref 3 provides a good summary of other 
laminar and turbulent flow foreed convection beet transfer correlations, These correlations may also be adapted for 
supercritical hydrogen cooling babes. 
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Hi. Property Uncertainty Analysis 


When conducting any experiment, die measurement* of all variable* have uncertainties associated with them, as 
do die values of material properties that ate obtained through reference databases. There uncertainties propagate 
throughout the data reduction equation and produce an overall uncertainty for the calculated result, For the case of 
our heat transfer analysis, the data reduction equations considered are Eqs. (2X (2>X (5X (EX and (8X with 





( 12 ) 


where the property values are determined at the reference conditions given hy Eqs, (4X (7) and (1 OX 

Property values for hydrogen were obtained from the NIST Thermodynamic and Transport Properties qfPure 
Fhtidt database 34 and the NIST Chemistry WehBoak,*' There sources are based on Ref. 22 for density and specific 
heat, and Ref. 23 for viscosity and thermal conductivity, They contain data for temperature up to 400 K, For 
conditions above 400 K, Ref 24 was used, which is an updated version of Ref 2$ to correct fer errors in certain 
property data. 

The uncertainties hsfed in Table 2 are compiled horn those supplied hy the referenced sources and are based on 
accuracies of experimental measurements and inaccuracies introduced in fitting the data, The uncertainties quoted 
are not meant to represent the critical region where the uncertainties could reach 10 to 20% or lugher. Indeed, an 
analytical function such as the Modified Benedict- Webb-Rubin Equation used by Ref, 22 cannot represent the 
proper behavior near the critical point, and is not considered to he valid in that region. For fee purpose of this 
analysis, fee uncertainty in the Mr correlation is not calculated near the critical region. 

To determine fee total uncertainty in M* resulting from fee uncertainties of each individual property variable, a 


general uncertainty analysis was performed, 
This was accomplished by using a Monte 
Carlo simulation program feat was created in 
C++. To perform fee Monte Carlo simulation, 
random errors ware selected for each 
variable, based on a uniform parent 
population wife a equal to one-half of fee 
assumed uncertainty for each variable at the 
particular temperature and pressure (as per 
Table 2) A simulated measurement was 
produced by summing fee drawn error wife 
the assumed “true” value of fee variable This 
was performed for each variable, and fee 
resulting Nu calculated at fee specific 
temperature and pressure This simulation 
was repeated 5,000 times at each temperature 
and pressure wife fee resulting uncertainty 
af-Yw determined as twice fee calculated 
sample standard deviation. 3 * The simulation 
was repeated fer each combination of 
temperature and pressure within fee range of 
interest. The uncertainties in experimentally 
measured variables ( T,P . etc.) were not 
considered here 



Figure 2: Flaw Chart of Analysis Program and Monte Carlo 
Simulation 
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Table 2: Summery of Property Data, Uncertainty Vataee and Sources Used la Analysis 



•Where ranges of uneertslntie* we provided, • Uim*t appronimetien I* used to represent die uneertsinty In the snslysls, 


The results of foe uncertainty analysis are shown in Figures 3 and 4- The blackened region on each plot is the 
near critical region where the uncertainties are not defined, In this region the resultant uncertainty ii believed to be 
highest, but no reasonable estimate exists to quantify it precisely. As ean be seen faun the figures, the McCarthy & 
Wolf and Taylor correlations, both based on bulk reference type, produce comparable uncertainties, The uncertainty 
is lowest at approximately 1,8% under the Uqwd-vepor line, 2% for J>< 3, and 5 to 1%, for 3 < TF< 12, The 
reduction in uncertainty at Tr> 12 to approximately 4% is a result of switching the source of property data in this 
region. The Hendricks correlation shows an uncertainly value of about 5,4% for 7V< 2, increasing to a maximum of 
7% at higher temperatures, The Hess & Kura and Miller correlation results reveal considerably higher uncertainties 
for these correlations at low pressures and temperatures and along foe liquid-vapor line. The uncertainty is highest at 
approximately 10%, and ia upwards of 7 to 8% to Pr ■ 2. For higher temperatures, the uncertainty is similar to that 
of the Hendricks correlation at 5 to 7%, 
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Figure 4. Uncertnlntlei Associated with Nu In Supercritical Region Sue to Uncertelntlei In Property Seta. 
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IV. Experimental Data Study 

To SIMM the possible effect* of the property uncertainty, it i« useful to compare experimentally determined 
heat transfer coefficients to these calculated using the heat transfer correlations. With that to mind, fee authors 
gathered published experimental heat transfer data for hydrogen in the supercritical regimes, and reduced that data 
using mere recent hydrogen property databases. 

The data used to this report originate from experiments conducted by Hendricks. Simoneau. and Friedman 
(NASA Lewis, J96SX* Hendricks, Graham, Hsu and Friedman (NASA Lewis, 1966), 18 and Aerojet-Genera) 
Corporation (1967>** The ranges of conditions tested to these experimental programs can be found to Table 1. 
These sources were selected because their datasets were readily available, and because together they represent a 
wide range of fluid property conditions, Combined, they represent a total of 2824 data points, 

figure 6 shows the reduced data in tetms of the deviation of ^ , where is determined using the various 

heat transfer correlations by the relation, 


K* 


(13) 


was calculated by (he respective investigators at each potaU^/^* *1 represents an exact agreement between 
predicted and experimental results, 

As can he seen to the figure, the McCarthy & Wolf and Taylor correlations (both hulk reference types) over- 
predict fee heat transfer in fee region close to the critical point. The Hendricks, Hess & Kunz, Miller and Schaet & 
Quenhneyer all under-prediet in this region. 


Sralfli Prill 



figure 8, Comparison of Predicted Heat Transfer Coefficient ef Cerrefetfeus hr Various 
Reduced Pressures. All plots are forRe t * 2e6, T t ■ 300 K , JDO.Olm , x/D m 40 , 
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Figure 6. Deviation of h^/h^ from Unity. 
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Tins tread can be seen in a theoretical comparison of the correlation* shows in figure 5 (b) for reduced pressure 
equal to 1 ,6, The spike in die prediction of McCarthy & Wolf and Taylor is therefore shown to not represent die true 
conditions w this region. 

Also from Figure 6, it can he seen that the McCarthy & Wolf and Taylor correlations under-predict in die region 
of the line of IV *1 far FV>4. As the pressure is increased, this range of underqpredtotioa extends fa larger values 
of Tr up to about Tr m 2, This can be seen in Figures 5 (e) and (d) where the trends for McCarthy & Wolf and 
Taylor can be seen fa diverge increasing from the other correlations at low temperatures and higher pressures. The 
Schaet & Quentmeyer appears to give die most accurate prediction close to the 7>*1 line for Pr <10, 

Plots of die experimental data points are shown graphically compared fa each correlation in the Appendix, The 
dashed lines on these plots represent die range of ±100% from the prediction of the correlation, From these plots, 
along with Figure ti, it is evident fast the deviation of the experimental results from those predicted by the 
correlations are in many eases considerably larger than can be explained by the effects of property uncertainty alone, 

V. Coadwdons 

1, The uncertainty in Nusselt correlations for supercritical hydrogen as a result of uncertainty in property data 
ranges from about 1,8 to 10%, 

2, The uncertainty on Nu is effected by the reference temperature type used, In the low pressure and low 
temperature range, toe bulk reference type correlations have the lowest uncertainty, 

3, The deviation of experimental results to those predicted by correlations cannot he explained solely by die 
uncertainty in property data 

4, The deviations are greatest near the erideal region and close to die transition line (7>*1) from liquid to 
supercritical gas for all pressures- The deviation along the transition line is important far cooling channel 
applications, and an improved correlation should be developed to better predict heat transfer in this region, 

5, The best overall fit to the experimental data eame from the Taylor (lq 8t) and Schaet & Qeuntineyer (Eq, 9) 
correlations, The intepal property type of Schaet & Quentmeyer is promising, and substituting it into other 
general Nusselt-type equations could lead to improved correlations 
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